Phase separation and dynamical arrest for particles interacting with mixed 
potentialsthe case of globular proteins revisited 
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We examine the applicability of the extended law of corresponding states (ELCS) to equilibrium 
and non equilibrium features of the state diagram of the globular protein lysozyme. We provide 
compelling evidence that the ELCS correctly reproduces the location of the binodal for different ionic 
strengths, but fails in describing the location of the arrest line. We subsequently use Mode Coupling 
Theory (MCT) to gain additional insight into the origin of these observations. We demonstrate that 
while the critical point and the connected binodal and spinodal are governed by the integral features 
of the interaction potential described by the normalized second virial coefficient, the arrest line is 
mainly determined by the attractive well depth or bond strength. This article is published in Soft 
Matter. The reference is: DOI: 10.1039/c0sm01175d 
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I. INTRODUCTION 

Gel formation in colloidal suspensions has continued to 
attract considerable attention during the last few years. 
Particular attention has been given to particles interact- 
ing via a short-ranged attractive potential, and the role 
of the key parameters volume fraction <j>, interparticle in- 
teraction strength Ua and range of the potential D in 
determining the state diagram has been addressed [1,2]. 
Depending upon the strength of the attraction, they can 
form amorphous solids from arbitrarily low volume frac- 
tion, </>, up to random close packing. 1 At intermediate val- 
ues of U a /kT and <f>, one encounters an intriguing inter- 
play between spinodal decomposition and gelation that 
results in the phenomenon of an arrested spinodal decom- 
position and the formation of a solid-like network [26]. 
While there are a number of experimental studies that 
report on the phenomena, there are relatively few quan- 
titative results on the actual location of the arrest line 
in the vicinity of the binodal and spinodal. In a recent 
study Lu et al.6 provided for example evidence that gela- 
tion for short range attractive particles indeed occurs as 
a consequence of an initial equilibrium liquidliquid phase 
separation. The authors claimed that all phase separat- 
ing samples gelled, i.e., that the gel line coincided with 
the binodal in the Baxter sticky sphere model, and that 
the concentration of the resulting colloidal gels would al- 
ways be around <f> <~ 0.55. This was interpreted as a 
result of the gel line intersecting the phase separation 
boundary at (f> ~ 0.55, with the dense phase undergoing 
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an attractive glass transition after phase separation. 

Moreover, these authors also suggested the existence 
of a universal state diagram when using the second virial 
coefficient B 2 , normalized to the second virial coefficient 
Bhs = 27T(7 3 /3 of hard spheres with the same diameter 
o , as an effective temperature B 2 = BijB^g. It had in 
fact already previously been suggested that the so-called 
extended law of corresponding states (ELCS) [7] should 
not only allow for a rescaling of thermodynamic proper- 
ties, but that it should also hold for dynamic properties 
such as gel or glass lines that describe the occurrence of 
dynamical arrest [2,8]. 

Another example, where a state diagram that com- 
bines liquid liquid phase separation and dynamical ar- 
rest has been reported, was based on the use of glob- 
ular proteins as model colloids [4,9]. Direct analogies 
between colloids interacting via a short range attractive 
potential and globular proteins have been used in the 
past quite successfully to understand the phase behavior 
of globular proteins and to rationalize some of the phe- 
nomcnological observations made during protein crystal- 
lization [1014]. In particular lysozyme has been exten- 
sively investigated in view of the potential but also lim- 
its of such a coarse-grained colloid approach to protein 
solutions [11,13,15,16]. 

While the existence of dynamical arrest in protein so- 
lutions has in fact already been reported previously, [13] 
it is only recently that we have been able to locate the 
glass line in the region below the coexistence curve for the 
globular protein lysozyme [4]. These findings arc sum- 
marized in the lysozyme state diagram in Fig. 1, where 
the temperature T is plotted as a function of the protein 
volume fraction (f>. A lysozyme solution above the co- 
existence curve quenched down to temperatures within 
region I undergoes complete macroscopic phase separa- 
tion after some time [4]. However, deeper quenches into 
region II lead to the formation of a gel due to the in- 
terplay between phase separation and dynamical arrest 
[4,9]. 
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Arrest is thought to occur during the early stage of 
spinodal decomposition when the volume fraction 02 of 
the dense phase intersects the dynamical arrest thresh- 
old (j>2,Giass given by the arrest line, which extends deep 
into the unstable region (right triangles in Fig. 1). At 
this point the dense phase forms an attractive glass and 
arrests, and phase separation thus gets pinned into a 
space spanning gel network with a characteristic length 
£ [4]. Additional support for this scenario was subse- 
quently obtained through investigations combining video 
microscopy and ultra small-angle light scattering experi- 
ments [9]. Moreover, this interpretation is also in agree- 
ment with the work of Dumetz et al.,[16] who looked at 
the phase behavior of several globular proteins and noted 
the existence of an arrested state. 




FIG. 1: State diagram of lysozyme suspended in an aque- 
ous Hepes buffer at pH of 7.8, with [NaCl] = 500mM. Shown 
are fluidcrystal coexistence curve (circles), binodal (squares), 
spinodal (triangles) and dynamical arrest line (right trian- 
gles), and the tie line separating region I from II. In region 
I we observe phase separation, in region II arrested spinodal 
decomposition, and region III corresponds to the attractive 
glass (see ref. 4 for details). 

If the invariance of the state diagram for colloids with 
short range attractions in the frame of a Noro-Frcnkcl 
ELCS would also hold for non-equilibrium properties 
such as gel and glass formation as proposed by a num- 
ber of authors, [6,8,17] this could serve as an ideal basis 
for comparing various systems and generating a generic 
phase diagram as a function of a single control parame- 
ter £>2- However, while there exist simulations that have 
addressed this problem, [8] there are to our knowledge 
no systematic experimental investigations available that 
have tested the applicability of the ELCS to glass or 
gel formation for colloidal suspensions. Moreover, while 
most theoretical, simulation and experimental studies 
have concentrated on purely attractive systems, any at- 
tempt to extend this approach to protein solutions and 
many colloidal systems of practical importance will have 
to consider the effects of amixed potential that also in- 
cludes a soft repulsion due to residual charges. While it 
has already been observed that the location of the bin- 
odal changes with ionic strength in a way that is at least 



qualitatively in agreement with the ELCS, [10] no infor- 
mation on the arrest line in the vicinity of the binodal has 
been reported so far, despite the fact that this is for ex- 
ample an essential ingredient of the zone picture for pro- 
tein crystallization. We have thus started a systematic 
characterization of the binodal and the arrest line un- 
der conditions where we increase the contributions from 
the screened Coulomb repulsion. We in particular con- 
centrate on the question whether the resulting liquidgas 
coexistence curves and the arrest lines can be superim- 
posed by using as a scaling parameter also if the elec- 
trostatic repulsion becomes more prominent and longer 
ranged. 



II. MATERIAL AND METHODS 

We used the globular protein hen egg white lysozyme 
(Fluka, L7651) in 20mM HEPES buffer at pH = 7.8, 
where lysozyme carries a net charge of +8e.l8 We varied 
the electrostatic repulsion created by the residual charges 
of lysozyme at pH=7.8 by adding different amounts of 
NaCl. Details about the sample preparation are given 
elsewhere [4]. The buffers used in this study contain 200, 
300, 400, and 500mM NaCl. Taking into account the 
counter-ions from the proteins, this corresponds to a de- 
crease of the Debye length l D from about 19% of the 
lysozyme diameter for 200mM to about 12% for 500mM 
NaCl for a lysozyme solution at c =20 mg mL _1 .19 For 
the determination of the state diagrams at different ionic 
strengths we again used a combination of visual observa- 
tion, optical microscopy, static light scattering, rheology 
and centrifugation as described in ref. 4. The volume 
fraction of lysozyme f was calculated from the protein 
concentrations c as measured by UV absorption spec- 
troscopy using <f> = c/p, where p = 1.351 g cm 3 is the 
protein density. 

The second virial coefficient for lysozyme solutions 
was measured using static light scattering experiments 
on concentration series of dilute solutions of lysozyme. 
The light scattering experiments were performed with 
a commercial goniometer system (ALV/DLS/SLS-5000F 
mono-mode fiber compact goniometer system with ALV- 
5000 fast correlator) at a wavelength of 514. 5nm and 
a fixed scattering angle of 90°, leading to a scattering 
vector q = 0,022 nm. Figure 2 summarizes the results 
from our determination of the second virial coefficient 
B2 function of temperature for four different val- 
ues of the ionic strength. Here we used the relation 
Kc/R = l/M + 2N A B 2 c/M 2 , where K is the contrast 
term, c is the concentration, M is the molecular weight, 
and Na is the Avogadros number, to determine B 2 from 
the measured Rayleigh ratios R. Plotted in Fig. 2 are the 
values of the second virial coefficients obtained and nor- 
malized by their value for hard sphere, B 2 S — 2tt<t 3 /3. 
Also shown are the values of the critical temperature 
T c for the different salt concentrations. It is interest- 
ing to note that for all salt concentrations investigated, 
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the Bf = B2/B? s values at T c are 5^ = 2.55 ± 0.12 




5 10 15 20 25 30 

T (°C) 

FIG. 2: Temperature dependence of the normalized second 
virial coefficient B2/B| /s obtained from light scattering ex- 
periments for different ionic strengths. Lines are guides to 
the eye. Vertical dotted lines highlight the value of the criti- 
cal temperature T c for a given value of the ionic strength. 



III. RESULTS 

We first determined the location of the coexistence 
curves as a function of salt concentration. The resulting 
binodals are shown in Fig. 3A.We find that the criti- 
cal volume fraction <fi c = 0.17 ± 0.02 is salt independent 
and consistent with values from literature [10,20]. The 
critical temperature T c , however, increases with increas- 
ing salt concentration, reflecting the fact that the attrac- 
tive part of the interaction potential which drives phase 
separation becomes more dominant as the electrostatic 
repulsion becomes increasingly screened. 

It has already been demonstrated in various studies of 
liquidgas coexistence and crystallization in protein solu- 
tions that the ionic strength dependence of the location 
of the liquidsolid and liquid liquid phase boundaries can 
be rationalized at least qualitatively using an integral 
quantity of the strength of the interaction potential as 
given by B 2 [11,13,16,2125]. We have thus also deter- 
mined the second virial coefficient B 2 as a function of 
temperature for all ionic strengths with static light scat- 
tering (see ESI for additional information). We found that 
the normalized B 2 at T c is B 2 = 2.55 ± 0.12, identical 
within experimental errors irrespective of ionic strength, 
in qualitative agreement with the idea that equilibrium 
properties are dictated by integral features of the inter- 
action potential.We can in fact rescale the corresponding 
phase boundaries for all values of the ionic strength us- 
ing £>2 as an effective temperature as shown in Fig. 3B. 
This provides us with a direct verification of the applica- 
bility of the ELCS for lysozyme as a model for particles 
interacting via amixed potential. 

In a next step we now focus on the dynamical arrest 
line as a function of ionic strength and temperature. We 
determined their location by a combination of rhcology, 



microscopy and centrifugation experiments as described 
in detail in ref. 4. The results are also summarized in Fig. 
4. In contrast to the strongly salt-dependent coexistence 
curves, the location of the dynamical arrest line seems to 
be independent of ionic strength, and all the arrest lines 
at different ionic strengths overlap within the experimen- 
tal error bars. This holds not only for the points within 
the unstable region that were determined by a combi- 
nation of a rapid temperature quench and subsequent 
centrifugation, but also for the intersection with the bin- 
odal (point A in Fig. 1) determined independently by 
rheology and microscopy. 
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FIG. 3: Lysozyme state diagram showing the metastable bin- 
odal (open symbols) and the dynamical arrest lines (full sym- 
bols) for different ionic strengths, (a) T vs. <f> presentation: 
the dynamical arrest lines for different values of the ionic 
strength all overlap. Lines are guides to the eye. (b) T e // 
vs. <p presentation: in agreement with the ELCS, the binodal 
lines all overlap. 

These experimental observations thus lead to the fol- 
lowing scenario: the coexistence curves are not sensi- 
tive to the details of the interaction potential and can 
be rescaled using an integral quantity the second virial 
cocfficicntthat characterizes the overall potential. How- 
ever, this drastically changes when it comes to dynamical 
arrest. Here it appears that the short-range part of the 
potential, i.e. the temperature-dependent attractive well, 
determines its location. 

We next turn to the idealized mode coupling theory 
(MCT) [27] in an attempt to gain some qualitative in- 
sight into the trends observed experimentally. MCT re- 
quires the static structure factor S(q) as input, which we 
obtain from PercusYevick (PY) theory [28] applied to an 
interaction potential of double Yukawa form, U (r) / kT = 
-Uo/kTlKie-^-^+K^-^-^/KKi+K^r] for sep- 
aration distances r > 1 and complemented by a hard core 
for r < 1. Here we chose to use a simple mixed potential 
that would allow us to capture the general trends upon 
a variation of the ionic strengths. The parameters of the 
screened Coulomb part of the interaction, K 2 > and z 2 , 
have been calculated from the experimental parameters 
using the linear superposition principle of the linearized 
PoissonBoltzmann equation. Their (weak) temperature 
dependence has been neglected. In addition, K\ < and 
Z\ = 50, so that the first term in the interaction po- 
tential models a short-range attraction and the contact 
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value provides a suitable temperature, viz. —kT/Uo, for 
a qualitative comparison with the experimental data of 
Fig. 3. However, we refrain from attempting to model 
the noncquilibrium transitions quantitatively because of 
the limitations of MCT for very short-range attractions. 8 
Fig. 4A shows the arrest line observed within MCT as 
a function of added monovalent salt. At the higher vol- 
ume fractions increasing the ionic strength has the effect 
of shifting the glass transitions to higher <p, in agreement 
with ref. 26. However, at sufficiently low temperatures 
and particle concentrations the glass lines merge and be- 
come essentially independent of the ionic strength. In 
this region the location of the glassy arrest transition 
produced by the theory is found to be governed by the 
contact value of the interaction potential. Indeed, at 
comparable f the location of the glass line is unaffected 
by changes in the ionic strength (see Fig. 4A). MCT 
thus captures the experimentally observed trends. This 
is further demonstrated in Fig. 4B, where the theoretical 
predictions are scaled by B\. Whereas PY calculations 
for the spinodals reveal that they superpose reasonably 
well under these conditions (data not shown), the glass 
lines from the theory do not scale with B\ , confirming 
that under the prevailing conditions it is the well depth 
or bond strength, rather than B\, that governs the struc- 
tural arrest. 
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FIG. 4: Attractive glass transition determined from MCT 
for [NaCl]=oo (pure attraction), 500 mM, and 200mM NaCl. 
Plotted is the inverse contact value of the interaction potential 
(a) and the reduced second virial coefficient (b) as a function 
of the particle volume fraction. Calculations have been made 
on the basis of a double Yukawa interaction, as described in 
ref. 26, some examples of which are shown in the inset for 
different ionic strengths. 

We are thus left with the discrepancy between the re- 
sults from ref. 8, where evidence was presented that 
the ELCS would also be valid for dynamical arrest, and 
our own findings. However, it is important to point 
out that the simulations performed in ref. 8 were per- 
formed for purely attractive potentials of very short range 
e = A/a < 0.01. On the other hand, lysozyme possesses 
a mixed potential, and has been modeled with a range 
for the attraction of approximately 0.04 < e < 0.34, de- 
pending upon whether an isotropic potential or a patchy 
description has been used [15]. It is clear thatMCT pro- 
vides qualitative agreement between experiment and the- 



ory only. However, despite its shortcomings, MCT should 
be able to reproduce dynamical arrest for particles inter- 
acting via a not too short-ranged potential reasonably 
well, [29] and indeed captures the trends observed (Fig. 
3 and 4). 



IV. CONCLUSION 

Additional simulations and experiments with colloidal 
model systems under comparable conditions (mixed po- 
tential, not too short-ranged attraction) are thus needed 
before we can generalize our findings. However, we be- 
lieve that our results are of considerable importance for 
current attempts to link proteins and colloids. A coarse- 
grained approach to the phase behavior of globular pro- 
teins based on a close analogy between colloids and pro- 
teins has frequently been used to gain insight into phe- 
nomena such as protein crystallization or protein con- 
densation diseases [10,13,14,30]. The assumption of a 
generalized state diagram has for example been used to 
develop new strategies in protein crystallization. How- 
ever, our study clearly demonstrates that there may be 
no sound basis for such a generic state diagram for glob- 
ular proteins that incorporates both equilibrium as well 
as non-equilibrium features. 

We may also speculate about the importance of our 
findings for the colloid community. For colloids with 
mixed interaction potentials we may observe a signifi- 
cant shift of the arrest line with respect to the binodal 
at different ionic strengths. For weakly screened electro- 
static repulsion this could result in a shift of the binodal 
below the arrest line and thus open up a direct path 
to gelation that does not compete with a simultaneous 
spinodal decomposition. This is a scenario that has been 
discussed previously in attempts to generalize our cur- 
rent understanding of gelation and glass formation, but 
no convincing experimental evidence has been provided 
so far.2Given the enormous structural and mechanical 
differences found in colloidal gels and glasses that have 
formed via different routes, this opens up new avenues to- 
wards the formation of solid-like complex fluids that may 
be extremely interesting for applications in materials or 
food science [31]. 
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